Association of polygenic scores for autism with volumetric MRI phenotypes in cerebellum and brainstem in adults

Previous research on autism spectrum disorders (ASD) have showed important volumetric alterations in the cerebellum and brainstem. Most of these studies are however limited to case-control studies with small clinical samples and including mainly children or adolescents. Herein, we aimed to explore the association between the cumulative genetic load (polygenic risk score, PRS) for ASD and volumetric alterations in the cerebellum and brainstem, as well as global brain tissue volumes of the brain among adults at the population level. We utilized the latest genome-wide association study of ASD by the Psychiatric Genetics Consortium (18,381 cases, 27,969 controls) and constructed the ASD PRS in an independent cohort, the UK Biobank. Regression analyses controlled for multiple comparisons with the false-discovery rate (FDR) at 5% were performed to investigate the association between ASD PRS and forty-four brain magnetic resonance imaging (MRI) phenotypes among ~ 31,000 participants. Primary analyses included sixteen MRI phenotypes: total volumes of the brain, cerebrospinal fluid (CSF), grey matter (GM), white matter (WM), GM of whole cerebellum, brainstem, and ten regions of the cerebellum (I_IV, V, VI, VIIb, VIIIa, VIIIb, IX, X, CrusI and CrusII). Secondary analyses included twenty-eight MRI phenotypes: the sub-regional volumes of cerebellum including the GM of the vermis and both left and right lobules of each cerebellar region. ASD PRS were significantly associated with the volumes of seven brain areas, whereby higher PRS were associated to reduced volumes of the whole brain, WM, brainstem, and cerebellar regions I-IV, IX, and X, and an increased volume of the CSF. Three sub-regional volumes including the left cerebellar lobule I-IV, cerebellar vermes VIIIb, and X were significantly and negatively associated with ASD PRS. The study highlights a substantial connection between susceptibility to ASD, its underlying genetic etiology, and neuroanatomical alterations of the adult brain. Supplementary Information The online version contains supplementary material available at 10.1186/s13229-024-00611-7.


Introduction
Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder characterized by difficulties in social communication, restricted interests, and repetitive behaviors [1].Many observation about neuroanatomical differences in ASD have been identified using magnetic resonance imaging (MRI) [2][3][4].In particular, the role of the cerebellum in social cognition and the pathophysiology of ASD has been well-established [5,6].ASD-related genes co-expressed at the cerebellum in the developing brain have been observed to contribute directly to the pathogenesis of ASD [7].Cerebellar grey matter (GM) volumes correlate with core ASD symptoms [8].Beyond the core symptoms, the various associated symptoms of ASD affiliate with the disruption in the brainstem [9,10].The relationship of ASD with MRI volumetric changes in the cerebellum and brainstem in a genetic context is still underexplored at the adult population level.
ASD is highly genetically heterogeneous and affects ~ 2% of children globally [11].The genetic liability for ASD influences structural and functional changes in brain anatomy [12,13], and neuroimaging genetics of ASD has been proposed as a potential novel biomarker for the early diagnosis of ASD [14].The significant genetic influence in ASD is evidenced by the heritability estimates ranging from 70 to 90% [11,15], and single nucleotide polymorphism (SNP) based heritability of 11% [16].Polygenic risk scores (PRS) allow taking into account the cumulative totals of susceptible SNPs that are associated with specific target traits, with each SNP assigned a weight based on its effect size in association with the trait determined by genome-wide association studies (GWAS) [17].The latest GWAS on ASD enables the construction of ASD PRS that have been effectively explored in association with neuroanatomical changes such as cortical thickness [18,19], surface area, gyrification [19], and grey matter concentration [20].
A recent meta-analysis on cerebellar volumetric changes in ASD reported a weak but significant (uncorrected) larger global cerebellar volume, and vermal lobules I-V, VI-VII to be most reported with changes in ASD [21].A population-based cohort addressing children reported higher levels of autistic traits that were inversely associated with cerebellar cortex volumes [22].Contrarily, decreases in the left cerebellum [23,24] and reduced cerebellar volume [25,26] were also reported.In recent studies, no significant differences in the cerebellum were reported exploring a European cohort (274 individuals with ASD and 219 control) [27] and, examination of cerebellar grey-matter volume in children with neurodevelopmental disorders revealed no significant differences in volume of the cerebellar lobules in ASD compared to controls [28].While the findings are inconclusive, the MRI volumetric studies of ASD are mostly based on children or adolescents [2][3][4]28], and brain changes associated with age have been widely observed in ASD [3,4].One of the consistent findings in ASD is an accelerated brain growth in early childhood and deceleration in late childhood and adolescence [29][30][31].
In the last ~ 15 years, there have been hardly any studies published investigating the cerebellar, or brainstem volumetric changes in association with ASD among adults [3,4,10].Recently, Yang et al. (2016) quantitatively estimated regional GM volume abnormalities in adults (18-60 years) with ASD using a meta-analytical approach.The authors found reduced GM volume in the cerebellum with a significant association of age, sex, and intelligence quotient (IQ) [32].In a longitudinal study by Lange et al. (2015), cerebral development at the whole-brain and the regional levels was compared between individuals with ASD and controls (ages 6-35 years), highlighting a significant difference in growth curves with atypical volume decline in ASD [33].
While previous studies have considerably advanced our understanding of the association between ASD and cerebellar and brainstem volumetric alterations, these studies included neuroimaging data primarily focused on children or adolescents, and were limited mostly to casecontrol studies.The genetic liability for ASD in association with neuroanatomical changes across the general adult population remains underexplored.Moreover, an emerging perspective postulates that ASD may manifest along a continuum, with a normal distribution of autistic traits evident in the general population, challenging the traditional binary approach of diagnosis in psychiatry research, by recognizing potential intermediate outcomes [34,35].This approach is highly suitable for bridging the research gap regarding the connection between the genetic predisposition to ASD and the volumetric measurement of the cerebellum and brainstem among adults in the general population.
In this study, we aimed to investigate the association between ASD PRS and the volumetric MRI phenotypes in the cerebellum and brainstem among adults at the population level.In addition to cerebellar and brainstem GM volume, we exploited the whole-brain MRI volumetric data to investigate the global volumes of the brain, cerebrospinal fluid (CSF), GM, and white matter (WM).Forty-four brain MRI phenotypes among ~ 31,000 participants from the UK Biobank (UKB) cohort were explored.First, we investigated sixteen Total Volumes, including total volumes of the Brain, CSF, GM, WM, GM of the whole Cerebellum, Brainstem, and ten regions of the cerebellum (I_IV, V, VI, VIIb, VIIIa, VIIIb, IX, X, CrusI and CrusII).Second, we investigated twentyeight Sub-regional Volumes of the cerebellum including the GM of the vermis and Left and Right lobules of each cerebellar region (Fig. 1).We used the largest available GWAS dataset for ASD (N = 18,381 cases and N = 27,969 controls) [16] to construct ASD PRS for each subject.The utilization of increased GWAS sample size for PRS estimation, coupled with the large UKB sample of MRI measurements, robustly enhanced statistical power for examining the association between volumetric MRI phenotypes and ASD PRS within the European population.

Participants
We employed baseline data obtained from the UKB project, one of the largest prospective population-based cohorts in Europe of > 500,000 subjects aged between 40 and 69 years recruited in 2006-2010 [36].Our analyses focused exclusively on participants of European ancestry, as determined by self-report and multidimensional scaling (Data-field 22006), aligning with the ancestry criteria of the ASD GWAS conducted on individuals of European descent [37].Exclusion criteria included withdrawal of consent, non-caucasian origin, failure to meet standard quality control criteria [38], sex mismatch between reported and genetic data (with sex as a covariate to control its influence), and outliers in terms of heterozygosity and genetic relatedness (further details provided below).We excluded participants with incomplete MRI data.Participants were also excluded from the present analysis if they reported having received a diagnosis of dementia, parkinson's disease, multiple sclerosis, epilepsy, stroke, head or neurological injury, trauma, or other chronic neurological problems (nervous system infection, brain abscess, encephalitis, demyelinating disease, cerebral aneurysm, cerebral palsy, brain haemorrhage, or brain cancer) [39,40].Finally, our study model considered ~ 31,000 participants.Detailed participant enlisting illustrated in Supplementary Fig. 1.Ethical approval for UKB data collection was granted by the North-West Multicentre Research Ethics Committee (11/NW/0382), and the use of the data (under approved UKB application 30172) at our department was further approved by the Regional Ethics Committee of Uppsala, Sweden.Informed consents were obtained from all participants included in the study (https://biobank.ctsu.ox.ac.uk/crystal/field.cgi?id=200).

Image acquisition and pre-processing
The present study made use of imaging-derived phenotypes generated by an image-processing pipeline developed and run on behalf of the UKB [41].The T1-weighted structural MRI acquisition was performed with a 3T Siemens Skyra scanner (software platform VD13) equipped with a standard Siemens 32-channel RF receive head coil [42].Global volumetric MRI Fig. 1 Schematic summary of the study.We investigated the association between ASD PRS and forty-four brain magnetic resonance imaging (MRI) phenotypes among ~ 31,000 participants in the UK Biobank for whom quality-controlled data were available (Details in Supplementary Fig. 1).Primary analyses included sixteen MRI phenotype -total volumes of the brain, CSF, GM, WM, GM of whole cerebellum, brainstem and ten regions of the cerebellum (I_IV, V, VI, VIIb, VIIIa, VIIIb, IX, X, CrusI and CrusII).Secondary analyses included twenty-eight MRI phenotypes -sub-regional volumes of the cerebellum including the GM of the vermis and left and right lobules of each cerebellar region.L, left hemisphere; R, right hemisphere phenotypes were obtained using the FAST segmentation tool [43].Regional brain volumes were further extracted using the Harvard-Oxford cortical and subcortical atlases (https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/Atlases), as well as the Diedrichsen cerebellar atlas (http://www.diedrichsenlab.org/imaging/propatlas.htm).All volumes are reported in mm 3 .Detailed information on MRI data acquisition and processing is available at: https://biobank.ctsu.ox.ac.uk/crystal/crystal/docs/brain_mri.pdf.

Genotyping and genetic quality control
We utilized imputed genotyping data sourced from the UKB [44] in the study.Quality control (QC) procedures adhered to recommended protocols for PRS construction using UKB data [44,45].The analyses were restricted to SNPs meeting certain criteria: a minor allele frequency (MAF) > 0.5, an imputation information score (INFO) > 0.8, and genotype missingness < 0.02.Participants meeting genetic data QC standards within the UKB (Data-field 22020) [45] were included in the dataset.This involved excluding individuals with missing rates > 0.02 on autosomes, displaying sex discrepancies, heterozygosity outliers, genetic relatedness up to third-degree relatives, or individuals not identified as 'White British' based on genetic grouping (Data-field 22006) [45].

Polygenic score construction
PRS computation were conducted for each participant in the UKB using the clumping and thresholding algorithms implemented in PRSice-2 (version 2.3.5)[46,47].The latest iPSYCH-Psychiatric Genomic Consortium (PGC) ASD GWAS summary statistics (November 2017), comprising 18,381 individuals with ASD and 27,969 individuals from the general population of European ancestry [16] was used as the base dataset for PRS construction.Following current guideline [48], we clumped SNPs in linkage disequilibrium (LD) using an LD-based r 2 ≥ 0.1 and a 250 kb genomic distance.PRS for seven P SNP -value thresholds (P = 1, 0.75, 0.5, 0.25, 0.1, 0.01, and 0.001) were generated, and selected to achieve a balanced signal-tonoise ratio due to the highly polygenic nature of ASD [16].The number of SNPs at each threshold is enlisted in Supplementary Tables 1, and a visual representation of these thresholds is provided in Supplementary Fig. 2.

Statistical analyses of PRS associations with MRI phenotypes
We conducted linear regression analyses using standardized PRS as the independent variable (Eq.1).The fortyfour MRI phenotypes of interest were sorted into two separate analyses (Primary and Secondary) to explore structurally from overall to deep and, facilitate elaborative discussion of the findings.In the primary analyses, MRI phenotypes of interest were sixteen Total Volumes, including total volumes of the Brain, CSF, GM, WM, grey matter of the whole Cerebellum, Brainstem, and ten regions of the cerebellum (I_IV, V, VI, VIIb, VIIIa, VIIIb, IX, X, CrusI and CrusII).In the secondary analyses, MRI phenotypes of interest were twenty-eight Sub-regional Volumes of the cerebellum including the grey matter of the left and right lobules of each cerebellar region (I_IV, V, VI, VIIb, VIIIa, VIIIb, IX, X, CrusI and CrusII) and eight vermis areas (VI, VIIb, VIIIa, VIIIb, IX, X, CrusI and Cru-sII).Covariates included age, age 2 , sex, genotyping batch, and the first 15 genetic principal components (Data-Field 22009).The age is the one recorded at the time of participant recruitment in 2006-2009, alongside the baseline genotyping data.The imaging data were obtained from the imaging visits from year 2014 and onwards [44].Additionally, imaging sites/scanners (Data-Field 54) to control for multi-site scanner subtle effect [49] and, the x, y, and z coordinates of the head position in the scanner (Data-Field 25756-25758) to control for static-field heterogeneity were included [50].Total volumes of the Brain, GM, WM, and CSF were analysed and normalised for head size (mm 3 ).For the brainstem and all cerebellum phenotype regression models, intracranial volume (a composite of total volumes of Brain, GM, WM, and CSF) was considered as an additional covariate.All covariates were standardized.P-value significance (α) for the regression analyses (Eq. 1) was corrected for multiple testing including both, the number of regions (sixteen in primary analysis and twenty-eight in secondary analysis) and the number of P SNP -value thresholds using False-Discovery Rate (FDR) [51] at 5%.Statistical analyses were performed using R (version 3.6.3).MRI illustrations were made using MRIcroGL (version 1.2.20220720).

Results
We investigated the association of ASD PRS and fortyfour volumetric phenotypes including sixteen 'Total Volumes' and twenty-eight 'Sub-regional Volumes' on ~ 31,000 subjects (Supplementary Fig. 1).Participants included in the analysis comprised ~ 48% female and ~ 52% male, aged 40-70 years with a mean age of ~ 55 years.

Total volumes
Seven primary MRI phenotypes were significantly associated with the PRS for ASD at one or more P SNP value thresholds (FDR = 5%) as illustrated in Fig. 2, and listed in Supplementary Table 2. Regions showing the highest variance explained by PRS are illustrated in Fig. 3.The proportion of variance explained by PRS varied between metrics in the range of 0.018% ≤ R 2 ≤ 0.034%.

Sub-regional volumes
Three secondary MRI phenotypes, including the left cerebellar lobule I-IV, vermis VIIIb, and vermis X, were significantly associated with the PRS for ASD at one or more P SNP value thresholds (FDR = 5%) (Fig. 4, Supplementary Table 3), Regions showing the highest variance explained by PRS are illustrated in Fig. 5.

Sensitivity analysis
We repeated the association analyses between ASD PRS and all forty-four MRI volumetric phenotypes (Total volumes and Sub-regional volumes), excluding the individuals (N = 25) with self-reported physician-posed ASD diagnosis (Data-field 20544) in UKB.The results remained consistent and statistically significant as the main analyses (Supplementary Table 4).

Discussion
This study is the first to investigate the linkage between the genetic propensity of ASD and forty-four volumetric brain MRI phenotypes, with a focus on dissecting the association between ASD PRS and volumetric shifts in the cerebellum and brainstem in adults.Additionally, we directly examined the correlation between the genetic predisposition for ASD and the global volumes of the brain, GM, WM, and CSF in an adult population.Beyond prior evidence on the neuroanatomical atypicalities of ASD, this study robustly highlights neuroimaging findings in relation to the genetic propensity of ASD, at the population level in an adult European cohort with substantial statistical power.
Our investigation using PRS constructed from the recently published GWAS for ASD indicate that adults genetically predisposed to ASD exhibit a significantly reduced total brain volume, WM volume, and brainstem volume, alongside with an increased CSF volume.Specifically, individuals with higher ASD PRS show a reduced volume in the cerebellar lobules I-IV, IX, and X, as well as vermes VIIIb and X.Our findings intriguingly highlight changes in the mature brain of adults at risk for ASD, while previous studies mostly investigated young developing brains.Most of these studies reported increased volumes among participants with ASD and are based on children or adolescents [2][3][4].However, volumetric changes with age among individuals with ASD have been widely discussed [3,4].An accelerated growth leading to an enlarged brain in children with ASD has been observed during early life, while a decrease in brain volume is found over the developmental period in which typically developing children experience slight increases [29,31].There is also consensus that the underlying pathology of brain development in ASD involves increased growth of both GM and WM during early childhood, followed by a rapid synaptic pruning process during early adolescence [52,53].
We show that adults with a higher genetic predisposition for ASD are more likely to have a reduced total brain volume.Volumetric analyses in adults with ASD reporting total brain volume changes are rare.To the best of our knowledge, our study is the first to report a total brain volume reduction in adults at risk for ASD.Additionally, we show that adults who are genetically predisposed for ASD show reduced WM volumes.There are mixed reports on total WM volume in ASD brains, indicating overall or region-specific increases that have been observed [4].However, a significant reduction of the mean volume of cortical WM was reported in a longitudinal study comparing participants with ASD and typically developing controls, which showed a developmental curve reflecting reduction with age [33].Our findings regarding a reduced total volume of brain and WM in adults genetically predisposed to ASD are in line with these findings.It can be hypothesized that individuals with ASD continue to show a higher reduction in brain volumes with age compared to individuals without ASD.
We found a substantial correlation between reduced brainstem volume and elevated ASD PRS.Our findings support the most commonly reported structural alteration of the brainstem in ASD, as described in a recent review on ASD neuroimaging studies [10].The review summarized mixed results in relation to the brainstem size in ASD.However, a reduction of the total volume or of at least one of the brainstem components was most commonly observed [10].The contribution of an abnormal brainstem neurotransmission to ASD symptoms becomes well evident when co-occurring features of ASD such as sensorimotor challenges or psychophysiological atypicalities are examined [9].Moreover, the age effect in ASD brainstem volumetric change is highlighted in the mentioned review, as indicated by the difference in outcome between adult and child populations.Additionally, a reduction in this difference with an increase in age is observed in longitudinal studies [10].
Our findings on increased CSF volume are in line with previous reports on adult brain changes in ASD [54,55].A recent review paper concluded as well that CSF volume is increased in ASD [25], as hitherto mostly observed in children or adolescents.The consistent evidence of increased CSF volume in ASD [25,54,55], has not yet been thoroughly explored in relation to ASD pathophysiology.However, in a recent review Shen (2018) discussed the potential links between ASD and CSF abnormalities and concluded that CSF could be a potential early biomarker for ASD and its subtypes.Infants with ASD, compared to typical infants, exhibit factors such as an imbalance in CSF production and drainage, slower CSF recycling in early life, and a higher ratio of CSF volume in infancy.These collectively contribute to a reduced ability to eliminate inflammatory metabolites and toxins, Fig. 5 Anatomical localization of polygenic risk effects for ASD on the secondary MRI phenotypes, Sub-regional Volumes.(A) Coronal cerebellum sections, (B) sagittal sections and render view (right) depicting the significant association between the cerebellar sub-regional volumes and ASD PRS.The color bar indicates the highest variance (R 2 ) explained by the ASD PRS for the secondary MRI phenotypes with significant associations eventually leading to increased levels of Aβ in neurons, sleep problems, and cognitive impairment involving the lymphatic system [56].
Our investigation on the cerebellum demonstrates that adults genetically predisposed to ASD anatomically present reduced cerebellar volume in regions including the lobules I-IV, specifically in the left hemisphere, IX, and X, as well as the vermes VIIIb and X.Previously, a decreased GM volume in the anterior lobe of the cerebellum (lobule I-V) and lobule IX has been reported in children with ASD [8].Our findings robustly show a similar, reduced volume of regions I-IV and IX among adults with a genetic predisposition for ASD.While the lobule I-V is involved in increased stereotypical and repetitive behavior, and lobule IX is correlated with poorer communication skills as measured by standard ASD scales [8], we posit that the cerebellar differences observed in children remain at the adult age.
While we observed that a reduced GM of the cerebellar region X is significantly associated with a genetic propensity for ASD, hardly any studies investigated this region in adults with ASD.Recent study including 58 individuals with ASD, and 34 typically developing individuals aged 8-30 years reported a reduced volume in the cerebellar region X (right lobule), which was inversely associated with ASD symptom severity [57].Furthermore, postmortem studies also reported dysplasia and hypoplasia of the cerebellar region X, indicating cell migration abnormalities [58].Lower Purkinje cell density in region X is also possibly involved in restricted eye movement behaviour in ASD [59].Our findings on the association between the cerebellar region X volume and ASD PRS support these findings.
While investigating the cerebellar volumes at the subregional level, we observed that adults with a genetic propensity for ASD have a reduced volume in the left cerebellar region I-IV, and vermes VIIIb and X.There are very few studies on adult ASD investigating the subregional level of the cerebellum.Overall, the cerebellar vermis volume has been consistently observed to be reduced in ASD [25].However, our investigation of the cerebellar vermis region among adults suggests that areas VIIIb, and X are significantly associated with ASD PRS.Previously, an increased volume of vermis VIIIb among children with ASD was observed [8].However, we found a decreased volume among adults.Thus, our results support a previous report that showed a decreased GM in vermis VIIIb occurring in early childhood and no effect with increasing age [60].Reduced cerebellar vermis volume is associated with ASD symptomatology, including impaired global development, cognition, expressive language, and gross and fine motor skills, as well as behavior problems [61].
We also observed the GM volume of the left I-IV cerebellar region to be significantly negatively associated with the ASD PRS.Interestingly, reduced GM in the right I-IV cerebellar region has been observed to correlate with higher severity of results in ASD diagnostic interviews and social interaction scores in children with ASD [8].Altered lateralized neurodevelopment of brain has been reported in ASD [62], which explains the observed left-right asymmetry.Further longitudinal neuroimaging studies are required to explain the age effect on variations in hemispheric patterns of volumetric changes.While we observed significant associations between ASD PRS and the left cerebellar regions I-IV, as well as vermes VIIIb and X, the total volume of the cerebellum was not significantly associated with ASD PRS.A plausible interpretation of such finding is that genetic factors implicated in ASD may specifically affect the development or size of certain cerebellar regions rather than the entire cerebellum uniformly.It is intriguing to speculate that these regions are particularly sensitive to the genetic variants associated with ASD.
Our findings robustly infer the understanding that age has a significant effect on the developmental trajectory of ASD [3,4] and support the theory of age-specific anatomic abnormalities in ASD [29,31].Our study suggest the age-specific deceleration to be predominant in adult life in ASD, resulting in reduced volume of the brain on the overall and specifically on the cerebellar regional level.Our findings of volume reduction on the cerebellar sub-regional level robustly support the previous consistent findings in post-mortem analyses of ASD individuals, reporting reduced Purkinje cell size and number regardless of age, sex, or cognitive ability [5].
Strengths of this study include the application of a powerful PRS approach to investigate the association between genetic predisposition for ASD and brain structure in the adult population.We utilized the recent-largest GWAS for ASD including multiple P SNP -value thresholds for the inclusion of risk variants in the PRS construction.One of the most methodologically diverse and largest MRI datasets to date was used to investigate the PRS associations with brain structures.The findings of the present study should, however, be interpreted with caution, as factors like an insult or socio-emotional behavioural damage during life [63], IQ and behavioral phenotypes [64] of participants have not been taken into account in the study, which may have contributed to the regional adult cerebellar volumetric differences.A limitation of the study is that the ASD PRS only explains a relatively small variance of 2.5% compared to a SNP-based heritability of 11% in ASD [16].Advancements in future GWAS with even larger sample size may potentiate our understanding of the associations observed in this study.This must also be acknowledged that, the findings in the study are based on the genetic predisposition for ASD and not based on ASD diagnoses.While our findings are in line with and support previous findings, future studies in similar direction on the ground of diagnosed ASD adult cases may provide more emphasized and elaborated understanding.Furthermore, the UK Biobank represents an aging cohort predominantly of European ancestry, which is characterized by higher average socioeconomic status and better health compared to the broader population [65].Thus, the generalizability of the presented findings may need further explorations in more diverse populations.
In summary, this was the first study to investigate associations between brain volumes focusing on the cerebellum and brainstem, and polygenic risk scores for ASD in the adult population with high statistical power, utilizing both genetic and imaging data.We show that certain global volumes such as the whole brain, WM, brainstem, CSF, and volumes in cerebellar lobules I-IV, IX, and X, as well as vermes VIIIb, and X are significantly predicted by ASD polygenic risk scores.These findings show that the cumulative genetic load for ASD is associated with considerable changes in the adult brain.ASD typically manifests early in life and entails significant neurodevelopmental symptoms, often accompanied by structural changes in the brain.This supports a possible explanation for our observations, suggesting that genetic susceptibility to ASD may lead to notable structural variations in the brain, even among individuals who do not exhibit the full spectrum of the disorder.The findings may robustly be confirmed by future translational MRI and histological studies involving human post-mortem data or animal models of genetic susceptibility for ASD.

Fig. 2
Fig. 2 Associations between polygenic risk scores for ASD and the primary MRI phenotypes, Total Volumes.Bar charts of variance explained by ASD PRS (R 2 , y-axis) constructed at each of seven probability thresholds (0.001 ≥ P SNP ≤1, x-axis) for each of the sixteen total volumes: total volumes of the Brain, CSF, GM, WM, grey matter of the whole Cerebellum, Brainstem and ten regions of cerebellum (I-IV, V, VI, VIIb, VIIIa, VIIIb, IX, X, CrusI and CrusII).Blue bars indicate negative associations and red bars positive associations; asterisks indicate P values for association after FDR correction: *P ≤ 0.05, **P ≤ 0.01, ***P ≤ 0.001.Polygenic risk scores for ASD were significantly negatively associated with the total volumes of the Brain, WM, Brainstem, and Cerebellar regions I-IV, IX, and X; and significantly positively associated with the CSF volume

Fig. 4
Fig. 4 Associations between polygenic risk scores for ASD and the secondary MRI phenotypes, Sub-regional Volumes.Bar charts of variance explained by ASD PRS (R 2 , y-axis) constructed at each of seven probability thresholds (0.001 ≥ P SNP ≤1, x-axis), for each of the twenty-eight cerebellar regional volumes: grey matter of the Left and Right lobules of each cerebellar region (I-IV, V, VI, VIIb, VIIIa, VIIIb, IX, X, CrusI and CrusII) and eight Vermis areas (VI, VIIb, VIIIa, VIIIb, IX, X, CrusI and CrusII).Blue bars indicate negative associations and red bars positive associations; asterisks indicate P values for association after FDR correction: *P ≤ 0.05, **P ≤ 0.01, ***P ≤ 0.001.Polygenic risk scores for ASD were significantly negatively associated with the volumes of the cerebellar regions I-IV (left), VIIIb (vermis) and X (vermis)